home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zhecon.z / zhecon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZHHHHEEEECCCCOOOONNNN((((3333FFFF))))                                                          ZZZZHHHHEEEECCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZHECON - estimate the reciprocal of the condition number of a complex
  10.      Hermitian matrix A using the factorization A = U*D*U**H or A = L*D*L**H
  11.      computed by ZHETRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE ZHECON( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, LDA, N
  19.  
  20.          DOUBLE         PRECISION ANORM, RCOND
  21.  
  22.          INTEGER        IPIV( * )
  23.  
  24.          COMPLEX*16     A( LDA, * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      ZHECON estimates the reciprocal of the condition number of a complex
  28.      Hermitian matrix A using the factorization A = U*D*U**H or A = L*D*L**H
  29.      computed by ZHETRF.
  30.  
  31.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  32.      condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
  33.  
  34.  
  35. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  36.      UPLO    (input) CHARACTER*1
  37.              Specifies whether the details of the factorization are stored as
  38.              an upper or lower triangular matrix.  = 'U':  Upper triangular,
  39.              form is A = U*D*U**H;
  40.              = 'L':  Lower triangular, form is A = L*D*L**H.
  41.  
  42.      N       (input) INTEGER
  43.              The order of the matrix A.  N >= 0.
  44.  
  45.      A       (input) COMPLEX*16 array, dimension (LDA,N)
  46.              The block diagonal matrix D and the multipliers used to obtain
  47.              the factor U or L as computed by ZHETRF.
  48.  
  49.      LDA     (input) INTEGER
  50.              The leading dimension of the array A.  LDA >= max(1,N).
  51.  
  52.      IPIV    (input) INTEGER array, dimension (N)
  53.              Details of the interchanges and the block structure of D as
  54.              determined by ZHETRF.
  55.  
  56.      ANORM   (input) DOUBLE PRECISION
  57.              The 1-norm of the original matrix A.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZHHHHEEEECCCCOOOONNNN((((3333FFFF))))                                                          ZZZZHHHHEEEECCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      RCOND   (output) DOUBLE PRECISION
  75.              The reciprocal of the condition number of the matrix A, computed
  76.              as RCOND = 1/(ANORM * AINVNM), where AINVNM is an estimate of the
  77.              1-norm of inv(A) computed in this routine.
  78.  
  79.      WORK    (workspace) COMPLEX*16 array, dimension (2*N)
  80.  
  81.      INFO    (output) INTEGER
  82.              = 0:  successful exit
  83.              < 0:  if INFO = -i, the i-th argument had an illegal value
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.